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ABSTRACT 

The South Pole Telescope (SPT) has nearly completed a 2500 deg^ survey of the southern sky in 
three frequency bands. Here we present the first public release of SPT maps and associated data 
products. We present arcminute-resolution maps at 150 GHz and 220 GHz of an approximately 
95 deg^ field centered at R.A. 82.7°, decl. —55°. The field was observed to a depth of approximately 
17 /iK-arcmin at 150 GHz and 41 /iK-arcmin at 220 GHz during the 2008 austral winter season. Two 
variations on map filtering and map projection are presented, one tailored for producing catalogs of 
galaxy clusters detected through their Sunyaev-Zel'dovich effect signature and one tailored for produc- 
ing catalogs of emissive sources. We describe the data processing pipeline, and we present instrument 
response functions, filter transfer functio ns, and map noise properties. All data products described 



in this paper are available for download at http://pole.uchicago.edu/public/data/maps/ra5h30dec-55 



and from the NASA Legacy Archive for Microwave Background Data Analysis server. This is the first 
step in the eventual release of data from the full 2500 deg^ SPT survey. 

Subject headings: cosmology: observations — surveys — cosmology: cosmic background radiation — 
methods: data analysis 
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1. INTRODUCTION 

1.1. Signals in the Millimeter-wave Sky 

Millimeter-wavelength (mm-wave) maps of the sky 
contain rich cosmological and astrophysical information. 
Away from the Galactic plane, the mm-wave sky is dom- 
inated by the cosmic microwave background (CMB) at 
large and intermediate angular size scales. The CMB 
features on scales of roughly five arcmin to many de- 
grees primarily arise from temperature fiuctuations at 
the surface of last scattering, and mm-wave measure- 
ments of these anisotropics have enabl ed powerful con- 
straints on cosmological niodels (e .g., 'Komatsu et al.| 
2Qll[ |Dunkley et al. 2QTT| [Ke isleTet al 2011). Interac- 
tions oT CMB photons with matter can induce secondary 
anisotropics through processes such as gravitational lens- 
ing and scattering, and measurements of these additional 
anisotropics are potentially powerful probes of cosmic 
structure. The CMB features on scales smaller than 
about five arcmin are dominated by secondary anisotropy 
due to the inverse Compton scattering of CMB pho- 
tons by fre e electrons, known as the Su nyaev-Zerdovich 
(SZ) effect ( Sunyaev fc Zerdovich|1972 ). On these small 
scales, emissive extragalactic sources also contribute sig- 
nificantly to the mm-wave sky. These sources are of con- 
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siderable astrophysical and cosmological interest them- 
selves, particularly the population of high-redshift, dusty, 
star- forming galaxies (DSFGs) that ma ke up the bulk 



of th e cosmic infrared background (GIB, [Lagache et aL 
|2005"). The current generation of mm- wave telescopes is 
just beginning to exploit the scientific potential of these 
small-scale signals. 

The SZ effect consists of two components: the kinetic 
SZ (kSZ) effect and the thermal SZ (tSZ) effect. The 
kSZ effect is due to Doppler shifting of CMB photons 
by the bulk velocity of electrons along the line of sight. 
The tSZ effect is due to the scattering of CMB pho- 
tons by hot, thermally distributed electrons, primarily 
in galaxy clusters. This interaction results in a spec- 
tral distortion of the CMB with a null at approximately 
220 GHz. At observing frequencies above this null, the 
tSZ effect produces an increment in measured CMB flux, 
while at frequencies below the null, it produces a decre- 
ment. The resulting tSZ features in the mm- wave sky 
can be used to detect and characterize massive galaxy 
clusters. Galaxy clusters trace the largest peaks in the 
matter density field of the universe, and their abundance 
as a function of mass and redshift is a sensitive probe 
of structure growth. Galaxy clusters selected in a fine- 
angular-scale tSZ survey provide a nearly mass-limited 
and nearly redshift-independent sample for constraining 
cosmological parameters such as the dark energy equa- 
tion of state parameter w and t he normalization of th e 
matter power spectrum as (e.g., Carlstrom et al.|[2QQ2 ). 
Such constraints are already being realized with just the 
first small fraction of data from large mm- w ave surveys 
dVanderlinde et al.||2QlQ[ jSehgal et al.||2QlQ ). 

Measurements of the tSZ power spectrum provide ad- 
ditional cosmological constraints, independent of those 
from the primary CMB and from c atalogs of individually 
detected clusters (jLueker et al. 2010 ; Dunkley et al. 2011] 
[Shirokoff et al.||20ll[ ). These measurements can also in- 
form models of the physical processes in gala xy clusters 
(e.g., |Shaw et al ||2010j [Battaglia et al.||2010[ ). Measure- 
ments oTTEeTS^ power spectrum are sensitive to, among 
other pro cesses, t he reionization history of the universe 
(e.g., Za hn et al.|20 05), and limits on kSZ from mm- wave 
measurem ents have already ruled out some reionization 
scenarios dMortonson fc Hu||2010| ). 

Two populations of extragalactic sources contribute 
most significantly to the mm-wave sky away from the 
Galactic plane: sources for which the fiux decreases or 
is roughly constant with frequency, consistent with syn- 
chrotron emission from active galactic nuclei ( AGN) ; and 
sources for which the fiux increases with frequency, con- 
sistent with thermal emission from DSFGs. Since the dis- 
covery by the SCUBA instrument (Holland et al. 1999) 
of a population of moderate-to-high-redshift DSFGs re- 
sponsible for a significant fraction of the total GIB emis- 
sion, DSFGs have been an active area of mm-wave and 
sub-millimeter (sub-mm) research. The recent discovery 
of a sub-population of str ongly lensed, h igher-redshift 
DSFGs ( Vieira et al.^2010i [Negrello et aL][2Q10} has fur- 
ther increased interest in these sources. Measurements 
of synchrotron-emitting sources in mm-wave bands have 
the potential t o constrain models of AGN physics (e.g.. 



The South Pole Telescope (SPT) is a 10-meter tele- 
scope designed to survey a large area of the sky at mm 
and sub-mm wavelengt hs with arcminute angular resolu- 
tion and low no ise (Ruhl et al.| |2004| [Padin et aL||2008 



Carlstrom et al.||2011| ). The current ^PT receiver is a 
three- band (95, loO, and 220 GHz) bolometer camera op- 
timized for studying the CMB and the tSZ effect. Since 
the SPT was commissioned in 2007, the majority of ob- 
serving time has been spent on a survey of 2500 deg^ 
designated the SPT-SZ survey. The final data set from 
the SPT-SZ survey will consist of maps for 19 contiguous 
subfields of 70 to 230 deg^, observed during the austral 
winter seasons of 2008 through 2011. The final depth 
for most of the survey will be approximately 42, 18, and 
85 /iK-arcmin^^ at 95, 150, and 220 GHz, with roughly 
200 deg^ (including the field discussed in this work) hav- 
ing deeper 220 GHz data. Following completion of the 
SPT-SZ survey, the receiver will be reconfigured for po- 
larization sensitivity and will image a subset of the SPT- 
SZ survey area to significantly lower noise levels. 

While three-frequency data from the complete 
2500 deg^ survey will eventually be released, the first 
public map release from the SPT-SZ survey consists of 
data taken during the 2008 season, when the SPT re- 
ceiver was primarily sensitive in the 150 and 220 GHz 
bands. This first release presents maps covering approx- 
imately 95 deg^ observed in those two bands. The area 
covered by these maps is referred to as the ra5h30dec- 
55 field, named for the J2000 coordinates of the approxi- 
mate field center. This was the first large field mapped to 
survey depth by the SPT and was centered on a 45 deg^ 
subregion optically surveyed by the Blanco Cosmology 
Survey24 (BCS, Desai et al., in prep.). The BCS data 
were collected before the SPT was deployed, in antici- 
pation of using the combination of optical and mm-wave 
data for joint galaxy cluster analyses. The center of the 
BCS region and, hence, the center of the ra5h30dec-55 
field, is R.A. 82.7°, decl. -55.0°. The ra5h30dec-55 
field has been studi ed in detail to extract power spec- 
trum measurements (Lueker et al."2010; iHall et al.|2010 ' 



Shirokoff et al. 2011; Keisler et al. 2011) and to produce 
catalogs of emissive sources (Vieira et al. 2010) an d SZ- 
selected galaxy clusters (^Staniszewski et aL ^ 2009 Van- 



derlinde et al. 
TEIi 



2010 



Williamson et al 
tHe" 



De Zotti et aL]|2010p 



1.2. The South Pole Telescope and Survey 



T MTr 

paper presents the 2008 SJrT maps of the 
ra5h30dec-55 field at 150 GHz and 220 GHz, discussing 
in detail the aspects of the instrument response and data 
processing that are relevant for interpreting and using the 
maps. We present maps with two variations in filtering 
and map projection, one designed for cluster-finding and 
one designed for the detection and characterization of 
emissive sources. It is not possible to summarize all rele- 
vant features of the maps in a simple set of data products 
without some loss of information, which limits the use of 
these maps for certain types of analysis. In particular, 
the data products presented here are not sufficient to pro- 
duce an extremely accurate noise model, nor to perform 
a cross-spectrum analysis for estimating the CMB power 
spectrum, nor to perform jackknife analyses to test for 

23 Throughout this work, map signal and noise amphtudes are 
expressed in units of K-CMB, expressing deviations from the av- 
erage measured intensity as equivalent temperature fluctuations in 
the CMB. 

http:/ /cosmology.illinois.edu/BCS 
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contamination (e.g., Stiirokoff et al. 2011 Keisler et al. 
2011). All of these would require all maps of individ- 
ual observations of the ra5h30dec-55 field, which is an 
order-of-magnitude larger set of data products than what 
we are currently releasing. Individual-observation maps 
may be included in future releases. 

Regardless of the intended use of the maps, we em- 
phasize the importance of understanding how the SPT 
instrument response, data processing, mapmaking, and 
noise properties affect the signals of interest. To that 
end, after describing the instrument and observing strat- 
egy in Sections 2]andj3) we focus the majority of the pa- 
per (Sections [4 and pTon a detailed discussion of these 
properties of the data. The maps themselves are pre- 
sented in Section |6) which also presents cross-checks of 
analyses using these maps compared to previously pub- 
lished SPT analyses. Section [7| describes the data prod- 
ucts available online. We provide an example calculation 
using these data products in Appendix [Aj 

2. INSTRUMENT 

The SPT is a 10-meter off-axis Gregorian telescope lo- 
cated at the National Science Foundation's Amundsen- 
Scott South Pole Station, where atmospheric conditions 
are among the be st in the world for mm and sub-mm 
observations (e.g., Radford 2011). The receiver images 
the sky using an array of transition-edge-sensor (TES) 
bolometers read out using frequency-multiplexed Super- 
conducting Quantum Interference Device (SQUID) am- 
plifiers. A description of t he instrument de sign and per- 
formance can be found in Carlstrom et al. (2011). Here, 
we summarize the aspects of the instrument design and 
performance that are most relevant for understanding the 
data products in this release. 

The detector array in the SPT-SZ receiver is made 
up of six wedge-shaped sub-arrays, each of which has 
140 detector pixels configured to observe in one of the 
95 GHz, 150 GHz, or 220 GHz observing bands.^^ These 
three observing bands have been selected to coincide with 
"windows" of high atmospheric transmission, and to op- 
timize discrimination of the tSZ spectral signature. The 
observing bands are defined on the low-frequency end 
by a circular waveguide coupled to each detector, and 
on the high-frequency end by low-pass metal- mesh filters 
mounted above each detector wedge (Ade et al. 2006). 
Detectors across a given wedge have similar bana-pass 
profiles, with an average bandwidth of 35 GHz for the 
150 GHz band and 44 GHz for the 220 GHz band in the 
2008 receiver. Details on the measured bands for 2008 



are given in Section 4.1 and discussed further in Bleem 
et al. (in prep.). 

Data from each individual bolometer channel are digi- 
tized at 1 kHz, then digitally low-pass filtered and down- 
sampled to 100 Hz before being written to disk. We 
refer to the resulting data stream as time-ordered data 
(TOD). The optical time response of each detector is 
approximately described by a single-pole low-pass filter 
with time constants varying between 10 and 30 ms for 
different detectors. Measurements of th e time response 
functions are presented in Section |4.2| The measured 
time response is deconvolved from the data during anal- 

Each detector wedge has 161 potential bolometer channels, of 
which 140 are read out. 



ysis, and additional anti-aliasing filters are applied when 
the data are binned to create maps, as described in Sec- 
tion [521 

Each detector's beam, defined as the response of the 
detector as a function of angle to a point source on the 
sky, is determined by the combination of a conical feed- 
horn abov e the detector and the optical design of the 
telescope (Padin et al. 2008). The main lobes of the 
beams are well described by two-dimensional Gaussians 
with average full-widths at half maximum (FWHM) of 
1.15 and 1.05 arcmin at 150 GHz and 220 GHz, respec- 
tively. Individual detector beam profiles vary, primarily 
depending on their placement in the focal plane. How- 
ever, a single average beam for each frequency is appro- 
priate for characterizing the effective angular response 
in the final maps. Beam measurements are presented in 
Section [44l 

Noise in the SPT-SZ data comes from four main 
sources: 1) noise due to statistical fluctuations in pho- 
ton arrival time, 2) noise intrinsic to the detectors, 3) 
noise from the readout system, and 4) brightness temper- 
ature fluctuations in the atmosphere, mostly due to in- 
homogeneous mixing of water vapor. The first two com- 
ponents are expected to contribute essentially "white" 
noise (i.e., equal amplitude at all temporal frequen- 
cies). Readout noise has a white noise component as 
well as a "1//" component, with power decreasing as 
temporal frequency increases. The spatial power spec- 
trum of atmospheric fluctuatio ns increases steeply with 
increasing spatial scale (e.g., Bussmann et al. 2005), 
leading to noise in the TOD that rises steeply at low 
temporal frequencies. The atmosphere dominates the 
SPT noise at frequencies below 1 Hz, while the pho- 
ton noise dominates at the higher temporal frequencies 
that correspond to the signal region for cluster and point 
source science. Typical single-detector noise-equivalent- 
temperatures (NETs) expressed in CMB units range 
from approximately 380 to 540 jiiK^/s at 150 GHz and 
640 to 850 flK^/s at 220 GHz for the data presented here. 
These NETs are estimated from the detector noise power 
spectra in a band corresponding to multipole i ~ 3000 in 
the scan direction, using the absolute temperature cali- 
bration described in Section [4.5.21 of this release. Given 
the typical numbers of detectors with good performance 
during 2008 observations, the total mapping speed was 
approximately 29 ^J.K^/s at 150 GHz and 74 /iKy^ at 
220 GHz for the 2008 season. 

3. OBSERVATIONS 

The primary observing mode for the SPT is scanning 
across the sky at constant elevation. Because the SPT is 
located within 1 km of the geographic South Pole, this 
corresponds almost exactly to scans at constant declina- 
tion. Complete observations of a field are assembled from 
many consecutive scans at stepped positions in eleva- 
tion. Throughout each roughly 36-hour cryogenic cycle, 
we perform multiple short calibration measurements in- 
terleaved with the field observations. These calibration 
measurements include observations of a small chopped 
signal from a non-aperture-filling thermal source, two- 
degree scans in elevation, and observations of the Galac- 
tic HII regions RCW38 and MAT5a (NGC3576), which 
are common calibration sources for mm-wave CMB ex- 



periments (Puchalla et al. 2002 Coble et al. 2003 Kuo 
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et al.||2QQ7 ). This set of regular calibrations allows us to 
characterize instrument response and monitor detector 
performance, as described in Section [4j 

Between February 13 and June 5 2008, 421 observa- 
tions were performed of the ~95 deg^ ra5h30dec-55 
field, each taking about two hours of observing time. 
Each complete observation comprised 176 constant- 
elevation scans across the field, with elevation offsets 
of 0.125 degrees between pairs of scans back and forth 
across the field. Twenty different initial starting eleva- 
tions were used for successive observations, at offsets of 
0.005 degrees. Variations in the starting position of suc- 
cessive observations enhance the uniformity of coverage 
in combined maps. Approximately half of the observa- 
tions of this field were performed with an azimuth scan- 
ning speed of 0.44 deg/s, with the remaining observations 
performed at 0.48 deg/s. 

4. CHARACTERIZATION OF INSTRUMENT RESPONSE 

4.1. Observing bands 

The SPT spectral bandpasses are measured using a 
beam-filling Fourier Transform Spectrometer (FTS), as 
described in Bleem et al. (in prep.). Transmission spec- 
tra were measured for ~50% of the detectors on each of 
the six detector wedges. For a given wedge, the detec- 
tor transmission spectra are highly uniform, with well- 
defined band edges at low and high frequencies that are 
set by a precision-machined circular waveguide and a 
common metal- mesh low-pass filter. For detectors in the 
same wedge, the band center has an rms variation of 
~1%. For each band, we construct an average response 
by weighting each detector's transmission spectrum by 
the inverse square of the detector's NET. Uncertainties in 
the final spectra in each band are dominated by the abso- 
lute frequency calibration of the FTS. We can verify this 
by comparing the measured vs. expected location of the 
low band edge (due to the circular waveguide cutoff), and 
we estimate this absolute frequency scale to be accurate 
to 0.3 GHz. The 150 and 220 GHz transmission spectra, 
averaged over all detectors in a given band, are available 
for download (see Section [t]), expressed as the response 
to a beam-filling, fiat-spectrum = constant) source, 
with the peak transmission normalized to unity. 

In Table [T] we give the band center and effective band- 
width for the 150 and 220 GHz bands. We have defined 
the band center to be 

I f{v)dv ' ^ ' 

where f(y) is the transmission spectrum averaged over 
all detectors in a given observing band, and the effective 
bandwidth is defined to be J f{u)du. 

To convert the SPT maps from CMB temperature 
units to intensity, one must consider both the spectrum 
of the source and the spectrum of the CMB. For a source 
with a spectrum = IoS{iy)^ this conversion factor is 

AT jAn{iy)S{iy)f{iy)du ' ^ ' 

where Art{h') is the telescope throughput, or etendue, 
and ^(i^, Tcmb) is the differential change in brightness 
of the CMB for a change in temperature. For beam- 
filling sources in a single-mode system such as the SPT- 
SZ receiver, An{u) = (? jv^. 



Spectra of astrophysical sources in mm-wave bands are 
typically approximated as power laws, such that l(y) = 
I^iylv^T. InTable[lj we give example conversion factors 
between CMB units and MJy/sr for beam-filling sources 
with OL values typical of some common mm-wave source 
families. The conversion factors are quoted for equal 
to the nominal band center, i.e., either 150 or 220 GHz. 
For comparison, we also quote the conversion factor that 
we would obtain if our bands were infinitely narrow and 
centered on the nominal band center. (The conversion 
factor in this case is simply ^(^, Tcmb) x 10^^, refiecting 
the definition of 1 MJy = 10-^° W m-^ Hz"^) 

To convert a measured temperature fiuctuation in a 
CMB map to an equivalent ther mal SZ Comptonization 



or Compton-?/ parameter (e.g., [Carlstrom et al. 2002 
one simply divides the measured by the mean CM 
temperature and a frequency-dependent thermal SZ fac- 
tor. For delta-function bands, this factor is equal to 



1 



4 (l + fe(^,Te)), (3) 



where x = hv/kB TcuB^ ^sz(^r^e) is ^ small relativis- 
tic correction (e.g., Nozawa et al. 2000j ), and Tg is the 
electron temperature of the cluster. For SPT bands and 
a beam-filling source, the effective band-averaged fsz is 
equal to 



(/sz) = 



!v-^'§{v,TcMB)f{v)dv 



(4) 



In Table [T] we give values of effective, band- averaged 
/sz for two values of Tg (0 and 8 keV) for both the real 
bands and the delta-function approximations. 

4.2. Detector time constants 

The temporal response function of SPT detectors over 
the signal band of interest for this work (< 10 Hz) can be 
described by a single-pole low-pass filter, with each de- 
tector characterized by a single time constant. The time 
constants are measured periodically using the chopped 
thermal calibrator. The time constants are estimated 
by fitting the amplitude and phase response of each de- 
tector using a sequence of chopper frequencies from 5 to 
10 Hz. Time constants are verified using fast scans of the 
detectors across bright astrophysical sources. Maps con- 
structed using only left-going scans can be subtracted 
from maps constructed with only right-going scans to 
verify that residual time constant errors do not con- 
tribute significant spurious signal in the combined maps. 
Exactly these tests were performed in CMB power spec- 
trum analyses of SPT data including the ra5h3Qdec-55 



field fcueker et al. 2010 Shirokoff et al. 2011 



et al-lzO llI), and no spurious signal was t'ouri3" 



Keisler 



le time 

constants do not change significantly over a season or 
with receiver temperature over the range of data used in 
the final maps, so for analysis purposes a single time con- 
stant parameter is associated with each detector for all 
observations during a given season. Time constants for 
the 2008 receiver configuration vary between about 10 
and 30 ms across the detector array, with median values 
of 19 ms and 17 ms for 150 GHz and 220 GHz detectors, 
respectively. 
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TABLE 1 
SPT Average Band Properties 
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real 


delta-fn. approx. 


real 


delta-fn. approx. 


Band center (GHz) 


153.4 


150 


219.8 


220 


Bandwidth (GHz) 


35.2 




43.7 




Conversion factors for power-law spectra 










Radio, a = -0.5 (MJy sr-^ R-^) 


396.3 


398.6 


476.7 


483.7 


Ray leigh- Jeans, a = 2 (MJy st~^ K~^) 


389.6 


398.6 


487.7 


483.7 


Dusty, a = 3.5 (MJy sr"! R-^) 


375.6 


398.6 


487.5 


483.7 


Effective /sz 










Te=0 


-0.923 


-0.954 


0.006 


0.038 


Te=8 keV 


-0.876 


-0.903 


-0.046 


-0.018 



4.3. Pointing Reconstruction and Astrometry 

The real-time pointing model used to control the tele- 
scope is initially calibrated using optical star cameras 
mounted on the telescope s t ructu re, as described in more 
detail in Carlstrom et al. (2011). The pointing recon- 
struction ort[ie"lnin-wave(^^ then calculated offline 
using daily measurements of Galactic HIT regions and 
information from thermal, linear displacement, and tilt 
sensors. 

During each cryogenic cycle, full observations are per- 
formed of the HIT regions RCW38 and MAT5a. For each 
of these observations, the response of each detector is fit 
to a scaled, translated version of a template image of the 
HIT region. After any modification of the focal plane or 
optical configuration (typically once a year) , a set of the 
RCW38 observations is used to measure each detector's 
pointing offset relative to the telescope boresight. The 
daily observations of both RCW38 and MAT5a through- 
out the observing season are used to constrain the point- 
ing model over time, and the RCW38 measurements con- 
tribute to estimation of relative calibrations, as described 



in Section 4.5.1 After the pointing model has been cor- 
rected using the HIT region and telescope sensor infor- 
mation, random errors in the pointing reconstruction of 
roughly 7^' rms remain from observation to observation. 
These random errors contribute to the width of the ef- 
fective beam in the final coadded maps, as described in 
the following section. 

The absolute astrometry of the final coadded maps was 
initially calibrated by comparing the SPT positions of 
a handful of sources to those sources' positions in the 
843 MHz Sydney Univer sity Molongolo Sky Survey cat- 
alog ( Mauch et al.||2QQ3 ). This calibration was expected 
to be accurate at the 10^' level. The recent publication 
of the t he Australia Telescop e 20 GHz Survey (AT20G) 
catalog (Murphy et al. 2Q1Q| ) provides an even more ac- 
curate astrometric calibration. The astrometry in the 
AT20G catalog is tied to VLBI calibrators and is accu- 
rate at the V level. Using 17 sources in the SPT 150 GHz 
data, we find very low scatter between SPT and AT20G 
positions. We do, however, see a small (< 10^') but statis- 
tically significant mean offset from the AT20G positions. 
We correct this offset by simply changing the definition 
of the field center for the ra5h30dec-55 maps from its 
nominal value of R.A. 82.70000°, decl. -55.00000° to 
R.A. 82.70247°, decl. -55.00076°. After this correction 
has been applied, the SPT positions of the 17 sources 
agree with the AT20G positions to better than 1" in the 
mean, with arcsecond-level scatter. An estimated un- 



certainty of 2" in each dimension accounts for statistical 
uncertainty in the SPT source detections and system- 
atic uncertainties due to potential offsets in the source 
centers at the ATCA and SPT observing bands and po- 
tential offsets between the SPT 150 GHz and 220 GHz 



maps. 



4.4. Beams 



A thorough understanding of the detector beams, or 
angular response functions, is critical for interpreting the 
signals in SPT maps. Sky signals are convolved with 
these functions in the process of observation, leading to 
the extended appearance of point sources in the maps 
and the suppression of small-scale power in the angular 
power spectrum. The structure of the SPT beams can be 
characterized as a Gaussian main lobe out to a radius of 
1 arcmin, near sidelobes at radii between 1 and 5 arcmin, 
and a diffuse, low-level sidelobe at radii between 5 and 
40 arcmin relative to the beam center. 

Dedicated observations of planets are valuable for mea- 
suring the far sidelobes, since planets are bright enough 
to adequately probe the tails of the response function. 
Planet observations are less useful for studying the main 
beam, however, because of the limited dynamic range 
of the detectors, and because the planets are extended 
sources. We use observations of the brightest quasar in 
each survey field to characterize the main beam shape. 
We reconstruct a two-dimensional profile of the beam 
by stitching together measurements of the "inner beam" 
within a 4-arcmin radius, and an "outer beam" covering 
radii from 4 arcmin to 40 arcmin. 

The outer-beam measurements are based on 7 dedi- 
cated observations of Venus performed in March of 2008, 
and 1 dedicated observation of Jupiter performed in Au- 
gust 2008. The observations consisted of a sequence of 
azimuth scans with 0.5 arcmin elevation steps between 
scans. Because the detectors are saturated when they 
observe the planets directly (and require several detec- 
tor time constants to recover), only data from the first 
half of each scan is used. For Jupiter observations, the 
impact of the Jovian satellites, which is very small to 
begin with, is mitigated by subtracting a template based 
on their known locations. The Venus and Jupiter data 
are filtered to remove atmospheric noise and CMB fluc- 
tuations, with the locations of the planets masked in the 
filtering. The average scan-synchronous signal, as mea- 
sured at distances larger than 40 arcmin from the planet, 
is subtracted from the maps. 

The inner beam shapes for this data release are 
measured using a bright quasar that appears in the 
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ra5h30dec-55 field itself and a bright quasar that ap- 
pears in the other field observed by SPT in 2008, the 
ra23h30dec-55 field. We have no evidence that the 
beam shape differs between these two fields. A small 
map is constructed around the brightest source in each 
field. The source is masked to a radius of 5 arcmin and fil- 
tering is applied to remove atmospheric noise and CMB. 
The residual CMB and noise in the central beam region 
leads to an approximately constant offset to the absolute 
response in this map. The inner beam profile based on 
the quasar map is then stitched together with the outer 
beam profile from Jupiter, with an iterative procedure 
that uses the Venus maps to determine the scaling be- 
tween the two components and to correct for the constant 
offset in the inner maps. 

This procedure leads to a composite two-dimensional 
beam profile that describes both the main beam and far- 
sidelobe response of the instrument. In addition, because 
the inner-beam measurement is based on the final coad- 
ded maps of the fields, it encapsulates the contribution 
of several-arcsecond random pointing variations to the 
effective angular response in the final maps. Pointing 
variations increase the effective beam width by approx- 
imately 3% and 5% at 150 GHz and 220 GHz, respec- 
tively. 

The inner and outer beam profiles for both bands are 
presented in Figure [l] The composite beam maps are 
available for download, as described in Section j?! 

In many applications, it is preferable to use a sim- 
plified approximation to the beam profiles rather than 
the full two-dimensional angular response functions. Az- 
imuthally averaged beam functions in map space and 
Fourier space are presented in Figure [2] Using the flat- 
sky approximation, we calculate the Fourier transform 
(FT) of the composite beam map, B{i^(j)£). From this, 
we compute the azimuthally averaged beam function. 




(5) 



There is a small bias (< 0.5% fractional error at ^ < 
10000) in this estimate of B£ due to residual map noise, 
and we remove this bias. 

The normalization of B{i) is somewhat arbitrary, in 
that it is degenerate with the absolute, CMB-power- 
spectrum-based calibration factor described in Section 
4.5.2 Our CMB-power-spectrum-based calibration uses 
the multipole range 650 < i < 1000, so we choose to nor- 
malize B{i) to 1 at ^ = 800 to minimize the correlation 
between beam uncertainty and calibration uncertainty. 

The uncertainty in our estimate of B{i) arises from sev- 
eral statistical and systematic effects, including residual 
atmospheric noise in the maps of Venus and Jupiter, and 
the weak dependence of Bi on the choice of radius used to 
stitch together the inner and outer beam maps. We con- 
sider seven sources of uncertainty in total. In Figure [2] we 
show the quadrature sum of these error estimates, which 
approximates the total uncertainty. The beam functions 
are uncertain at the few percent level, and this uncer- 
tainty increases mildly with increasing multipole num- 
ber. 

4.5. Calibration 



Averaging measurements from many detectors within 
a given observation requires correcting for their rela- 
tive gains. To average multiple single-observation maps 
together, we also need to correct for inter-observation 
detector gain variations and changes in atmospheric 
opacity that affect the calibration of the whole array. 
Section 4.5.1 describes our procedure for estimating 
detector-to-detector and day-to-day relative calibrations. 
When all of these relative calibration factors are applied, 
a single absolute calibration number relates the ampli- 
tudes in the final map to physical units. Section |4.5.2| 
describes our procedure for estimating the final absolute 
calibration for each map. The physical units chosen for 
our maps are K-CMB, expressing deviations from the av- 
erage measured intensity as equivalent temperature fluc- 
tuations in the CMB. 

4.5.1. Relative Gain Calibrations 

Two types of calibration observation are combined to 
account for detector-to-detector relative gains, temporal 
gain variations, and atmospheric opacity variations: 1) 
daily observations of the HII region RCW38, and 2) ob- 
servations of the chopped thermal calibration source that 
illuminates the focal plane from the center of the tele- 
scope secondary mirror. The thermal calibration source 
illuminates the focal plane through a small hole in the 
secondary mirror; hence, the effective source shape seen 
by the detectors is very different from that of sources on 
the sky, and varies over the focal plane. 

The thermal calibration source is observed many times 
per day, and we use these observations to correct for vari- 
ations in response to sky signal across the array (flat- 
fielding) and to correct for any gain drifts. The first 
step in the relative calibration pipeline is to relate each 
detector's response to the calibration source to its re- 
sponse to an astrophysical source. To do this, we as- 
sign an effective temperature to the calibration source 
for each detector, based on the season average of ratio 
between the response to the calibration source and the 
response to RCW38. The relative calibration for each de- 
tector over a single ra5h30dec-55 observation is then 
based on the response to the calibration source observa- 
tion nearest that field observation. To account for any 
day-to-day drifts in the calibration source filament tem- 
perature or illumination pattern, we correct the single- 
observation relative calibration number for each detec- 
tor using the wedge- averaged difference between that 
observation's calibrator response and the season aver- 
age. Similarly, to account for any changes in atmospheric 
opacity, we correct each detector's calibration using the 
wedge- averaged difference between response to the near- 
est RCW38 observation response and the season average. 
The relative calibration factors applied to a given detec- 
tor are typically stable to within 2% over the season. 

4.5.2. Absolute Calibration 

The HII regions RCW38 and MAT5a can be used as 
absolute calibrators, but their irregular shapes and the 
uncertainties in their absolute fluxes at mm-wavelengths 
limit the precision of the calibration. An RCW38-based 
absolute calibration has been use d for some previous 



work (Staniszewski et al. 2009[ [Williamson et al. 2011^ 
Story et al. 20lT] ), but here we adopt a refined absolute 



calibration tased on comparisons of CMB angular power 
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Fig. 1. — Two dimensional average beam functions for the ra5h30dec-55 field. Panels (a) and (b) show the 150 GHz and 220 GHz 
two-dimensional beam profiles on a linear scale, emphasizing the structure of the "inner beam." Panels (c) and (d) show the beam functions 
on a logarithmic scale out to much greater radii, emphasizing the structure of the "outer beam," and showing the stitching of the two beam 
estimates at a radius of 4 arcmin. Note that in panels (c) and (d) the absolute value has been taken in order to display the image on a 
logarithmic scale, but this visually exaggerates the appearance of the noise in the innermost 4 arcmin of the beam. 

spectra produced from SPT and WMAP data. The two 
approaches give calibrations that are consistent within 
their uncertainties. 

To calibrate the coadded maps to absolute temperature 
units, we estimate the angular power spectrum of the 



CMB at 150 GHz using the technique described in Keisler 
et al. ( 2Q11[) . The maps used in this estimation have been 
preliminarily calibrated using the RCW38-based calibra- 
tion. The resulting spectrum is then com pared to the 



well-calibrated WMAP7 power spectrum ( jLarson et al. 



2011 ) over the multipole range 650 < i < 1000, and the 
final calibration to apply to the 150 GHz SPT maps is de- 
termined by requiring the inverse- variance- weighted ratio 
of the SPT and WMAP7 power spectra to be unity over 
this range. To reduce the contribution of SPT noise and 
sample variance to the calibration uncertainty, we use all 
2008 data to derive the calibration for the ra5h30dec- 
55 field. We have looked for evidence of calibration dif- 
ferences between the 2008 data from the ra5h30dec-55 
field and the rest of that season's data by checking the 
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Fig. 2. — Azimuthally averaged beam profiles (left), Fourier-space beam functions (center), and fractional uncertainties (right). 



regular observations of MAT5a, and we constrain this dif- 
ference to be smaher than 1% in temperature. A cross- 
power spectrum is constructed using the 150 GHz and 
220 GHz maps, correcting for differences in beams and 
filtering between the two bands. The ratio of this cross- 
power spectrum to the power spectrum at 150 GHz for 
multipoles dominated by primordial CMB signal gives 
an estimate of the relative calibration offset between the 
two bands, which is used to transfer the 150 GHz abso- 
lute calibration to 220 GHz. 

We estimate that the uncertainty in our 150 GHz ab- 
solute temperature calibration is 3.1% using this calibra- 
tion method. The uncertainty in the 220 GHz absolute 
calibration is estimated to be 6.9%. Note that since the 
220 GHz calibration derives from the 150 GHz calibra- 
tion, the two uncertainties are correlated with a correla- 
tion coefficient of approximately 0.5. 

5. DATA PROCESSING AND MAP-MAKING 

Raw time-ordered data are processed into maps in two 
stages. In the first ("pre-processing") stage, the raw data 
from a single observation of the field are calibrated, data 
selection cuts are applied, and initial filtering and in- 
strument characterization are performed. The same pre- 
processing pipeline is used for all SPT survey data, re- 
gardless of the intended final use of the maps, although 
there have been some minor modifications in the data 
selection criteria and pre-processing algorithms over the 
history of SPT data analysis. 

In the second ( "map-making" ) stage, additional filter- 
ing is performed on the pre-processed TOD and the data 
are binned into single-observation maps used for final 
coadds. The filtering applied in the map-making stage 
is typically customized for different analyses, as is the 
choice of map projection. As described in detail below, 
we present two specific combinations of filtering and map 
projection in this release. 

We characterize the effect of the filtering and projec- 
tion choices by describing the Fourier-space properties 
of the maps. We refer to features in the Fourier plane 
using the corresponding angular wavenumbers in radi- 
ans k. The maps are oriented such that at the center 
of the map, the x-direction corresponds to R.A. and the 
^/-direction corresponds to declination. 

5.1. Data Selection and Pre- Processing 



The first steps in the pre-processing pipeline are to 
characterize detector performance, assemble measured 
response characteristics, and calculate calibration factors 
and weights. 

Summary statistics are first calculated to characterize 
the receiver setup and the TOD for each detector. Each 
detector's TOD for an entire 2-hour observation is used 
to estimate the power spectral density (PSD) for that 
detector. A fit to the noise PSD is then used to identify 
line-like features that deviate from a 1// plus white noise 
profile. Additional noise statistics are calculated for the 
SQUIDs. 

Next, a set of responsivity statistics is assembled for 
every detector. These include the amplitude of the de- 
tector's response to the chopped thermal calibrator, the 
amplitude of its response during a calibration observa- 
tion of RCW38, and the amplitude of its response to a 
2-degree elevation scan performed before each field obser- 
vation. These response parameters are used to calculate 
daily gain calibratio ns fo r each detector, as previously 
described in Section |4.5[ They are also used to define 
data selection cuts. 

Preliminary cuts verify that each bolometer's bias volt- 
age and readout configuration settings are within the 
nominal ranges and that the TOD values for each channel 
remain within the dynamic range of the digitizer. Once 
these initial cuts have been applied, bad bolometers are 
rejected with progressively more and more stringent as- 
sessments of responsivity and noise. 

The first such cuts enforce a minimum signal-to-noise 
response to the chopped calibration source, and a mini- 
mum response to the routine elevation scans. Bolometers 
are also rejected if their response to the elevation scans 
does not fit the expected modulation of the atmosphere. 
If a detector's response to the chopped calibration source 
or the short elevation scans is more than three standard 
deviations away from the median response for detectors 
of the same band, it is also rejected. Detectors are re- 
jected if their PSDs exhibit wide line-like features or too 
many lines. After all of these cuts are applied, a final pair 
of cuts rejects bolometers with calibration constants or 
noise weights that deviate by more than a factor of three 
from the median for each band. The median numbers of 
detectors that pass cuts for t he r a5h3Qdec-55 field ob- 
servations (quoted in Section 5.6) reflect the bolometers 
selected through this process. 

The bolometer data are parsed into scans, defined as 
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temporally contiguous periods of constant-velocity az- 
imuth scanning in either direction across the field. Cut- 
ting periods of time when the telescope is accelerating 
(at the ends of each scan) removes about 5% of the total 
TOD for an observation. All data for a given scan are cut 
if the receiver temperature exceeded the nominal operat- 
ing value at any time during the scan, if the telescope fol- 
lowing error (absolute value of the commanded position 
minus the position recorded by the encoders) exceeded 
20" at any time during the scan, or if there were data 
acquisition problems leading to bolometer data or point- 
ing data drop-outs during the scan. Typically about 5% 
of scans are cut for these reasons. Data for an individual 
bolometer in an individual scan are also rejected if that 
bolometer shows evidence of step-function features in its 
TOD (which can sometimes occur due to "fiux jumps" 
in the SQUIDs), if the SQUID or bolometer noise for 
that scan is excessive, or if any of the digitized SQUID 
data approach the limits of the dynamic range. A spike- 
finding algorithm identifies cosmic-ray-like events in the 
TODs for individual detectors. If there are fewer than 
five in a given scan, and the spike features are relatively 
small, we remove them and interpolate over the gaps; 
otherwise, we fiag and ignore affected TODs for the du- 
ration of the scan. Typically, around 5% of otherwise 
well-performing bolometers are cut from each scan for 
one of these reasons. 

Because a varying subset of detectors will sometimes 
show sensitivity to the receiver's pulse-tube cooler (a 
phenomenon also seen in ot her TES bolom eter systems 
with pulse-tube coolers, e.g.. Dicker et al.|[2009), we ap- 



time constants for every individual detector (see Sec- 
tion [42). A cut-off frequency of 25 Hz is used for the 



ply a notch filter to remove a small amount of bandwidth 
from all data during the pre-processing stage. Combin- 
ing all data that pass selection cuts for a given observa- 
tion, we identify the fundamental frequency of the pulse- 
tube cooler and notch-filter a conservative 0.007 Hz of 
bandwidth around this frequency as well as any strong 
harmonics. For the 2008 ra5h30dec-55 field observa- 
tions, the pulse tube operating frequency was 1.62 Hz. 
If all harmonics of this frequency were always cut, this 
would represent an absolute maximum of 0.4% of the to- 
tal bandwidth; the actual bandwidth cut is smaller since 
typically only the first few harmonics are cut. We have 
verified that this notch filter has a negligible effect on fur- 
ther analysis, and we neglect it in further analysis steps. 

After all pre-processing, the offline pointing model is 
used to calculate corrected pointing positions to asso- 
ciate with every time sample, and the pre-processed data 
are written into an intermediate data format for further 
analysis. 

5.2. Filtering 

After the pre-processing stage of data analysis, the 
TOD associated with a given field observation is fur- 
ther filtered on a scan-by-scan basis. The processing of 
each scan 1) deconvolves the detector temporal response 
functions, 2) removes high-frequency temporal noise that 
translates to spatial scales smaller than the pixel resolu- 
tion in the maps, and 3) removes atmospheric noise while 
minimizing the filtering of signal power. 

The time-constant deconvolution and low-pass anti- 
aliasing filters are applied in a single Fourier-domain 
operation on each scan for each detector. The time- 
constant deconvolution takes into account the measured 



low-pass filter applied to the TOD, chosen to limit noise 
on spatial scales smaller than the 0.25-arcmin map pix- 
els without suppressing power on the spatial scales of the 
SPT beams. 

Prior to averaging the data into maps, additional time- 
domain filtering is applied to remove atmospheric noise. 
The atmospheric contamination resides primarily at low 
temporal frequencies, so a high-pass filter is an effective 
way to remove the most contaminated data. However, 
the application of a Fourier- domain high-pass filter will 
also distort the appearance of bright point sources in the 
final maps, introducing "ringing" patterns. This effect 
can be avoided by filtering in the time domain via fitting 
to slowly varying template functions and masking the 
locations of known bright point sources during the fit. 
Point-source masking of this sort has been performed for 
one of the sets of maps in th is re lease, but not for the 
other, as explained in Section [53] 

In all variations of our analysis, we fit out a mean and 
slope as well as higher-order polynomials from each scan. 
In some past work we have subtracted Legendre polyno- 
mials and in others a series of sines and cosines (Fourier 
modes). In either case the result is a time-domain high- 
pass filter that allows for point-source masking as needed. 
For the maps in this release, we have fit Fourier modes 
up to a temporal frequency that corresponds to a spa- 
tial high-pass cutoff of k = 300 in the scan direction, or 
kx = 300 (because the scan direction corresponds to the 
X direction in our maps). 

Atmospheric noise is highly correlated across the de- 
tector array. By taking this correlation into account, we 
can filter additional atmospheric noise without affecting 
small-scale sky signal. For each TOD sample, we sub- 
tract the mean value across every detector in a given 
wedge. This filter can also be performed with the loca- 
tions of bright point sources masked, and this masking 
is performed in one set of maps in this release (as with 
the masking in the time-series filters). 

5.3. Filtering Variations for Maps in this Release 

This paper presents two sets of maps, one tailored to 
cluster analysis and one tailored to point-source analy- 
sis. For cluster analysis, the point sources in the maps 
are viewed as a foreground. The locations of known point 
sources are masked in the filtering for these maps in or- 
der to prevent artifacts from ringing, which could affect 
cluster extraction. The brightest 201 sources are masked 
in the filtering, with the top 9 masked to a radius of 5 
arcmin and the rest masked to 2 arcmin. The source list 
used for constructing the mask was created by merging 
together all sources detected at greater than 5<j signif- 
icance in either 150 GHz or 220 GHz ma ps, from the 
publi shed point-source catalog for this field ( Vieira et al. 
2010), and adding any sources that lie outside the field 
boundaries used for detection in that work but that are 
strongly detected in an analysis using all map pixels. 

When the point sources themselves are the object of 
study, it is preferable to filter them in the same way as 

The majority of extragalactic emissive sources appear point- 
like in the arcminute-resolution SPT maps, and we use "point 
source" and "emissive source" interchangeably in this work. 
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the other parts of the map, so that the effect of filtering 
on signal can be completely characterized. For the set 
of maps tailored to point source analysis, the brightest 
point sources are not masked in the filtering. We also 
choose different map projections for the two sets of maps, 
as described below. All other parameters of the data 
processing and filtering are equivalent in both sets of 
maps. 



5.4. Map Projections 

After all data selection and filtering has been per- 
formed, we use the corrected pointing information to av- 
erage and bin the data into pixels in a two-dimensional 
map. For small areas of the sky, analysis of the maps 
is greatly simplified by using a flat-sky approximation. 
This allows two-dimensional FTs to be used (rather 
than spherical harmonic transforms) for analyzing sig- 
nal power, noise, and instrument response as a function 
of angular scale. However, any choice of map projec- 
tion will lead to some distortions of information in the 
map. The character of these distortions varies depending 
on the selected projection, meaning that some analyses 
may be easier to perform on maps with one projection 
while others may be easier with another. 

In past work we have employed several different pro- 
jection schemes, of which two have been selected for pre- 
sentation here. Past SPT cluster-finding a nalysis has 
used the Sanson-Flamsteed projection (e.g., "Calabretta 
Greisen 20 02), which projects constant-elevation scans 
into pixel rows in the final map. This projection choice 
simplifies the characterization and treatment of signal 
filtering, because the TOD processing effectively oper- 
ates on map rows. The disadvantage of this projection is 
that it is not distance-preserving, and sky features near 
the corners of the maps will therefore have slightly dis- 
torted shapes (for example, a 1-arcmin circle will appear 
to have an ellipticity of e ~ 0.12 near the corner of the 
ra5h30dec-55 map). We have selected this projection 
for the maps tailored to cluster-finding. 

The second projection employed here, for the maps 
tailored to point-source analysis, is the oblique Lambert 
equal- area azimuthal projection. For the typical size, 
shape, and center location of SPT maps, this sky pro- 
ject ion2reseryesboth distances and areas to high accu- 
racy ( Snyder|1987 ), which means that the effective beam 
shape will not vary significantly for different locations 
in the map. This is particularly advantageous for re- 
constructing point sour ce amplit udes using variants of 
the CLEAN algorithm ( |HQgbom| [l974). The disadvan- 
tage of this projection is that the effective filtering has a 
strong position dependence, since the angle between the 
scan direction and the map pixel rows changes across the 
map. For constant-declination scans, this angle, a can 
be expressed as: 



and 
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a = tan"^(A/5), 

A = 7 sin ^0 sin(0o — 0) x 

(sin 00 cos ^ + sin ^ cos cos(^o — ^)) 
+ cos6>o sin((/)o - 0), 

5 = 7 sin ^0 sin^ (^o — ^) sin + cos(0o — 0) , 



(6) 



Here equals the R.A. in radians, while = 7r/2 — 
declination, also in radians. Unsubscripted variables rep- 
resent the pixel location and the subscript "0" denotes 

the field cent er. 

In Sections [6^ and [631 we will discuss the impact of the 
choice of projection on the signal and noise properties in 
these maps. 

5.5. Weights 

The weights used to coadd data from many detectors 
into a single-observation map are calculated by averaging 
the calibrated PSD for each detector between 1 and 3 Hz. 
For the telescope scan speeds used in observations of the 
ra5h30dec-55 field, this frequency range corresponds 
to a multipole range of roughly 1500 < £ < 4500 in the 
scan direction, a reasonable overlap with the scales of 
interest for source and cluster detection and for high-^ 
power spectrum analyses. 

When we average the data from individual detectors 
to obtain a measurement of the map value in a partic- 
ular pixel, we calculate the total weight for that pixel 
by combining the weights from each individual detector. 
These pixel weights are used to combine the multiple 
single-observation maps into a final coadded map (see 
Section 5.6). The primary source of noise variation be- 
tween single-observation maps is the weather, which is 
accurately reflected in the calibrated 1-3 Hz noise esti- 
mates for the individual detectors in all but the poorest- 
weather days. 

The weights for the final map are used to define the 
uniform-coverage region. We define the uniform-coverage 
region as the area of the resulting map at 150 GHz for 
which the weight in a map pixel exceeds 95% of the me- 
dian for all map pixels, resulting in an approximate area 
of 95 deg^. 

5.6. Coadding Single- observation Maps 

If all single-observation maps of the ra5h30dec-55 
field in a given observing band were identical up to noise 
variations captured accurately by the pixel weights de- 
scribed in the previous section, then the final maps in 
each band would simply be weighted averages over all 
the single-observation maps. In reality, some fraction of 
single-observation maps need to be excluded from the fi- 
nal coadd based on anomalous behavior in the telescope, 
the weather, or the detectors. 

Of the 421 total observations of the ra5h30dec-55 
field, we reject 55 from consideration for final maps be- 
cause they were performed with atypical gain settings for 
the detectors or had pointing control problems leading to 
atypical coverage of the field. We reject an additional 21 
observations that were incomplete, due to the need to 
cycle the cryogenic system. 

We examine the weights and the noise rms in the cen- 
tral five- by- five- degree region of the individual observa- 
tion maps and make further cuts based on the weight and 
noise rms statistics. We have found that two particular 
cuts are necessary to optimize the final map noise and 
avoid biased signal values in the final map. First, we 
cut any observations with anomalously high weights, a 
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Fig. 3. — Map of the ra5h30dec-55 field at 150 GHz, using point-source masking in the data filtering and employing the Sanson- 
Flamsteed projection. This is the filter and projection scheme tailored for performing cluster analysis. For display purposes, the map has 
been smoothed with a 1-arcmin FWHM Gaussian, and high-noise regions near the boundary of the map have been masked. 



condition which implies that the recorded detector noise 
rms is lower than we can reasonably expect, possibly be- 
cause changes in loading have significantly altered the 
detectors' operating point. We also cut individual maps 
in which the product of the median weight and the noise 
rms squared is anomalously high. This can occur if the 
total map rms on all scales does not track the 1-3 Hz 
noise on which the weights are based, possibly due to 
anomalously poor weather. The weights- and noise-rms- 
based cuts eliminate an additional 20 individual maps at 
150 GHz and 29 individual maps at 220 GHz. 

All observations of the ra5h30dec-55 field that sur- 
vive the data selection cuts were taken with the sun be- 
low the horizon. The moon was above the horizon for 
52% of these observations but was never closer than 80 
degrees from the field center. Cross-checks pe rformed 
for power spectrum analyses using these data ([Lueker 



et ai.|2010| |Shirokoff e t al.j2011[[Keisler et al.|20lT] ) have 
shown no evidence of contammation from moon pick-up. 
The observations of the ra5h30dec-55 field sample the 
full range of azimuth directions. Cross-checks performed 



for power spectrum analyses using the se data ([Lueker 
et al.|2010[^Shirokoff et al'. 2011 ; Keisler et al.|201ip have 



shown no evidence of contamination from ground-based 
signals on spatial scales that are measured in these maps. 

After all of these checks and cuts, the final 150 GHz 
maps presented here are built from 321 individual ob- 
servations totaling 621 hours, with a median number of 
304 individual detectors contributing data during each 
observation. For the 220 GHz maps, 313 observations 
totaling 605 hours are included, with a median number 
of 166 detectors contributing data. 

6. MAPS 

6.1. Overview 

SPT maps are constructed by a simple inverse-noise- 
weighted averaging of observations of a given map pixel. 
Noisy modes are filtered out of the timestream rather 
than de-weighted in the mapmaking, and the resulting 
map is not an unbiased estimate of the sky signal. An 
understanding of both the noise and the effect of filtering 
on signal is essential for interpreting the maps. 
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Fig. 4. — Map of the ra5h30dec-55 field at 220 GHz, using point-source masking in the data filtering and employing the Sanson- 
Flamsteed projection. This is the filter and projection scheme tailored for performing cluster analysis. For display purposes, the map has 
been smoothed with a 1-arcmin FWHM Gaussian, and high-noise regions near the boundary of the map have been masked. 



Ideally, the data processing and map projection would 
be optimized separately for different science goals. In 
order to present a relatively simple set of data products, 
we have selected two variations on filtering and projec- 
tion choices. The first set of maps, created with cluster- 
finding in mind, is produced using point-source masking 
in the TOD filtering, using the Sanson-Flamsteed projec- 
tion. The second set of maps, created with point-source 
characterization in mind, is produced without masking 
bright sources in the filtering, using the oblique Lambert 
azimuthal equal-area projection. 

The 150 GHz and 220 GHz maps with these process- 
ing and projection choices are presented in Figures [3] 
through [6j For display, each map has been smoothed 
by convolving it with a 1-arcmin FWHM Gaussian to 
suppress noise on scales smaller than the approximate 
size of the beam, and high-noise regions near the bound- 
ary of the map have been masked. Figure [7] shows 
one-dimensional (azimuthally averaged in £ space) sig- 
nal+noise and noise PSDs for each map (computed from 
the point-source-masked maps in the Sanson-Flamsteed 



projection). These plots are included to illustrate: 1) 
where in spatial frequency the maps are signal- or noise- 
dominated; 2) where the noise is white and where it has 
a "red" spectrum due to atmosphere; and 3) what sig- 
nals are contributing at different spatial frequencies. In 
the following sections we discuss estimates of the filtering 
and noise properties of the maps. 

6.2. Filter Transfer Functions 

The beam and the timestream filtering are the two re- 
sponse functions that have the most impact on the prop- 
erties of sky signals recorded in SPT maps. While these 
response functions are often combined, we separate them 
here because they vary in different ways with the choice 
of map projection. The beam functions were described 
in Section [T4I In this section, we discuss how the time- 
domain filtering of the data from each scan combines to 
affect the properties of signals in the two-dimensional 
maps. 

The data filtering, as discussed in Section [5^ involves 
low-pass and high-pass filters with a wedge-average sub- 
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Fig. 5. — Map of the ra5h30dec-55 field at 150 GHz, with point sources unmasked during data filtering, using the oblique Lambert 
equal-area azimuthal projection. This is the filter and projection scheme tailored for performing point-source analysis. For display purposes, 
the map has been smoothed with a 1-arcmin FWHM Gaussian, and high-noise regions near the boundary of the map have been masked. 

ferent as a point source and a simulated CMB sky. For 
the maps presented here, we have estimated the transfer 
functions using simulations of a Gaussian input signal 
with an FWHM of 0.75 arcmin, chosen to probe the spa- 
tial scales of interest in these maps. 

Figure [S] shows the estimated two-dimensional fil- 
ter transfer function for the maps tailored to cluster- 
finding (filtered with point sources masked and presented 
with the Sanson-Flamsteed projection) In the Sanson- 
Flamsteed projection, the telescope scan direction is 
nearly equivalent to the x direction in the map. This 
means that the low-pass and high-pass time-domain fil- 
ters translate effectively into low-pass and high-pass spa- 
tial filters along the x direction. The wedge-average fil- 
tering acts as a roughly isotropic high-pass filter. In ad- 
dition, it produces localized decrements in the transfer 
functions at the spatial scales corresponding to the sep- 
aration distances between individual detectors. Sky sig- 
nal modes with wavelengths at the detection separation 
scale and aligned in angle with the pixel configuration 
will be almost perfectly removed by the wedge-average 



traction at each time sample. These TOD filtering oper- 
ations combine with the coverage and projection to result 
in filtering of signal on certain spatial scales in the maps. 
We characterize the effect of filtering on sky signals by 
estimating the two dimensional filter transfer functions. 
These are the Fourier-domain functions representing the 
relative suppression of signal power as a function of an- 
gular scale in the x and y dimensions of the map. 

The transfer functions are typically estimated by sim- 
ulating observations of a known signal using the recon- 
structed pointing for the real data, and passing the simu- 
lated observations through the full data analysis pipeline. 
Dividing the two-dimensional discrete FT of the simu- 
lated map by the FT of the input signal yields an esti- 
mate of the filter transfer function. Assuming that the 
effect of filtering is linear, the estimated transfer func- 
tions should be independent of the input signal used in 
the simulations. In practice, the estimated transfer func- 
tion depends very slightly on the input signal used, but 
the difference averaged over all spatial modes of interest 
is less than 1%, even for input signals as drastically dif- 
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Fig. 6. — Map of the ra5h30dec-55 field at 220 GHz, with point sources unmasked during data filtering, using the oblique Lambert 
equal-area azimuthal projection. This is the filter and projection scheme tailored for performing point-source analysis. For display purposes, 
the map has been smoothed with a 1-arcmin FWHM Gaussian, and high-noise regions near the boundary of the map have been masked. 



subtraction. The arrangement of these features reflects 
the hexagonal symmetry of the detector array. The angle 
between the symmetry axis of the features and = 
reflects the fact that the focal plane is slightly rotated 
with respect to the vertical axis of the telescope (to in- 
crease uniformity of map coverage along the elevation 
direction). 

The width of the reconstructed map varies as a func- 
tion of declination in the projection, so there is a slight 
change in the effective spatial cut-off scales for the high- 
pass and low-pass filters at the top and bottom of the 
map. Neglecting these small effects, the filtering is es- 
sentially uniform across the maps. 

Filter transfer functions for the maps produced with 
the oblique Lambert equal-area azimuthal projection are 
more complicated. In this projection, the telescope scans 
are at an angle to the x-direction in the map, and this 
angle varies with position. This means that the two- 
dimensional transfer function rotates as a function of 
map position (the angle of rotation is equal to the an- 
gle between map rows and R.A., which is computed as a 



function of map position in Equation |6|. To address this 
rotation as well as more subtle changes with position, 
we have separately estimated the transfer functions with 
the input Gaussian signal placed in nine different loca- 
tions in the map (the changes in transfer function across 
the map are sufficiently slow and regular that nine loca- 
tions easily sample the full behavior). Figure |9] compares 
the estimated transfer function for the central portion 
of these maps to the estimated transfer function in the 
lower right corner of the maps. The corner transfer func- 
tion shows slight changes in the effective spatial scale of 
the low-pass filtering, and it exhibits additional diagonal 
features with amplitudes of 5-10%. These features arise 
because the filtered power associated with a given bright 
source is spread diagonally across multiple rows of the 
map. In other words, they are a result of the interaction 
between the filtering and the sampling function associ- 
ated with the map pixelization. For most applications, 
these features can be neglected, as long as the rotation 
is taken into account. The transfer functions for all nine 
sub-regions of the map are available online. 
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Fig. 7. — One-dimensional (azimuthally averaged in t space) signal+noise and noise PSDs for each observing frequency. Left Panel: 
signal+noise and noise-only PSDs from the raw map, uncorrected for beam and filtering effects. Right Panel: as in left panel, but with 
the beam and filter transfer function divided out. Azimuthal averages are calculated using rough noise weighting, with noisier modes at 
low kx (or ix) receiving less weight and modes below kx = <tx = 400 ignored entirely. These demonstrate that the 150 GHz map is signal- 
dominated at nearly all spatial frequencies out to = 10000, but that the 220 GHz map has significant noise contributions, particularly at 
t > 5000. The essentially white character of the map noise is evident above i = 2000 in the raw map PSDs. The signal+noise PSDs with 
the beam and transfer function divided out show that the damping tail of the primary CMB (steeply falling with i) dominates the signal 
below i '2^ 2500, while the contribution from point sources (fiat with i) dominates above i ~ 2500. 



6.3. Map Noise 

6.3.1. Map Noise Estimation and Description of Noise PSD 

Features 



Each individual two-hour observation of the field yields 
a map with nearly identical sky signal contributions but 
independent atmospheric and instrumental noise. We es- 
timate the two-dimensional noise PSDs for the final coad- 
ded maps using the jackknife noise estima tion technique 



(Sayers et al. 



2009 



2009 ) , which is used 



Halverson et al. 

and described in many previous 

Staniszewski et aL]|2009 Vanderlinde et aL]|2010| Vieira 
et al.||2010|). Here, we describe those functions, noting 



publications (e.g.. 



that we refer to the square root of the noise power spec- 
tral density function as the PSD. 

Figure [lO] shows the averaged two-dimensional noise 
PSDs estimated for the maps constructed with the 
Sanson-Flamsteed projection. For most of the signal re- 
gion of interest for cluster-finding and point-source char- 
acterization, the noise is essentially white. Because noise 
power in the TOD translates to slightly different spatial 
scales at the top and bottom of the map, there is an effec- 
tive noise gradient of about 10%, scaling with the square 
root of the cosine of the declination in the map. The 
PSDs sho wn in Figure [lO] and the depth values quoted 
in Section 6.3.2 are averages over all declinations in the 



map. 

The most obvious non-white features in the PSDs are 
the concentration of high noise at low k (near the ori- 
gin), the increase in noise toward low and the noise 
cutoff at kx values below ~ 300. The increasing noise at 
low kx is due to 1// noise that is uncorrelated between 
detector channels. The cutoff at low kx is due to the 
high-pa ss fi lter, just as in the filter transfer functions in 
Section [6^ High noise near the origin is due to bright- 
ness fiuctuations in the atmosphere. These fiuctuations 
are spatially isotropic and have a very "red" spectrum. 



with most of the power on large spatial scales. Depend- 
ing on the scan speed of the instrument, in any given 
observation the atmosphere will either be imaged as if 
it were a stationary sky signal — in which case the atmo- 
sphere will appear as isotropic noise at low /c, just as it 
does in the composite PSDs shown here — or it will blow 
past the array in the direction of the wind — in which case 
it will appear elongated in /c-space along the wind direc- 
tion in a single-observation PSD. Even in the latter case, 
changes in wind direction over the course of many obser- 
vations will tend to make the atmospheric contribution 
essentially isotropic. 

The visible decrements in the PSDs in Figure [lO] cor- 
respond to the features previously described in the two- 
dimensional transfer functions. Residual atmospheric 
noise contributes excess noise at large spatial scales, but 
the majority of the atmospheric and instrumental noise is 
removed by the high-pass filter. A small amount of noise 
power "leaks" into the central high-pass-filtered region of 
the Fourier plane because of the subsequent application 
of the wedge-average filtering. The areas of excess noise 
at high ky /low kx are explained as "upmixed" large-scale 
noise power. This upmixing occurs because the sampling 
and filtering of large-scale atmospheric features is slightly 
different in the TOD of each individual detector. 



Figure [TT] shows the averaged two-dimensional noise 
PSDs estimated for the center and lower right corner of 
the maps with the oblique Lambert equal-area azimuthal 
projection. The properties are similar to those for the 
Sanson-Flamsteed projection, with the added rotation 
for positions away from the center of the map (the rea- 
son for this rotation and the value of the rotation angle 
are identical to those discussed for the filter transfer func- 



tions in Section 6.2). In these PSDs, the effect of residual 



high-frequency readout-related line features is visible as 
vertical strips at high kx- These are less visible in the 
full-map average PSDs presented in Figure IlO) but the 
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(a) Filter transfer function for 150 GHz 



(b) Filter transfer function for 220 GHz 



Fig. 8. — Estimated two-dimensional transfer functions for the 150 GHz and 220 GHz maps employing point-source masking and the 
Sanson-Flamsteed projection. These are averages over the entire map. The dark stripe along kx = is due to the effective scan-direction 
high-pass filter from the polynomial and Fourier-mode subtraction, the dark spot centered at ky = kx = comes from the isotropic high- 
pass filter of the wedge mean subtraction, and the hexagonally spaced dark spots are due to the sensitiv ity o f th e we dge mean subtraction 
to modes with wavelengths and angles corresponding to the detector array configuration. (See Sections |5.2| and |6.2| for details on filtering 
and associated transfer function features.) 



features are present. These features have a neghgible 
effect on analyses performed with the maps. 

As with the transfer functions, PSDs have been esti- 
mated for nine sub-regions of the obhque Lambert equal- 
area azimuthal projection maps, and these data products 
are available to download. 

6.3.2. Map Noise Accounting 

The noise PSDs discussed in Section [6.3. II can be used 
to estimate the effective map depths, but it is important 
to highlight that these numbers depend on the absolute 
calibration and on the normalization of the beam func- 
tions. In previous SPT publications, map depth has been 
quoted under the approximation that the SPT beams 
are Gaussian, with B{i) = 1 dit i = 0. Using the mea- 
sured SPT beam presented here, with B{i) normalized 
to unity at ^ = 800, the quoted map depth numbers 
are slightly different, but, as we will show, the two sets 
of beam/calibration/depth numbers are consistent with 
each other. Put another way, the fundamental quan- 
tity of interest for science is the noise after the beam 
and transfer function have been deconvolved, and the 
deconvolved noise we measure here is consistent with the 
Gaussian-beam estimates in previous papers, as shown 
below. 

Quadrature- averaging an annulus of the PSDs between 
i = 4000 and £ = 5000, we estimate a depth of 14.8 /iK- 
arcmin at 150 GHz and 36.1 /iK-arcmin at 220 GHz. For 
this estimate, as in Figure [7| azimuthal averages are cal- 
culated using rough noise weighting, with noisier modes 
at low kx (or i^) receiving less weight and modes below 
kx = = 400 ignored entirely. This estimate of map 
depth must be interpreted using the beams reported in 



this release, because the calibration applied to the maps 
depends on the beam function shape and normalization. 

For simple calculations of signal sensitivity, it is com- 
mon to use a Gaussian approximation to the beam. If we 
approximate the beams as Gaussians with FWHM 1.15 
arcmin and 1.05 arcmin for the 150 GHz and 220 GHz 
maps, respectively, we also need to modify the effective 
calibration of the maps. Within the i range used for es- 
timating the map depths, the average ratio of the true 
beam function to the Gaussian beam function can be 
used to adjust the depth figures. Over the £ range used 
in this noise estimate, that ratio is 0.86 at 150 GHz 
and 0.87 at 220 GHz. The depths appropriate for a 
Gaussian approximation to the beam are thus approx- 
imately 17 /iK-arcmin at 150 GHz and 41 /iK-arcmin 
at 220 GHz, as quoted in the abstract. These numbers 
are consistent with the quoted values of 18 /iK-arcmin 
and 40 /iK -arcmin from previous S PT publica tions (e.g., 
yanderlinde et al.|[2QTQ| [Shirokoff et~al][20TT] ) . 

I'he 150 GHz beam-deconvolved noise at i = 3000 (see 
Figure [7| right panel) is 19.7 /iK-arcmin. This is equal to 
the beam-deconvolved noise amplitude at this multipole 
for a 1.15-arcmin FWHM Gaussian beam and 18 /iK- 
arcmin white noise. The 220 GHz beam-deconvolved 
noise ai i = 3000 (see Figure [tI right panel) is 46.9 /iK- 
arcmin. This is equal to the beam-deconvolved noise at 
this multipole for a 1.05-arcmin FWHM Gaussian beam 
and 43 /iK-arcmin white noise. 

In Section |2j we reported mapping speeds for the full 
detector array at 150 and 220 GHz, using the absolute 
temperature calibration described in this release. These 
numbers can be combined with the total integration time 
on the ra5h30dec-55 field reported in Section [3] to ob- 
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Fig. 9. — Estimated two-dimensional filter transfer functions for the 150 GHz and 220 GHz maps employing no point-source masking 
and the oblique Lambert equal-area azimuthal projection. The "center" filter transfer functions represent the innermost one-ninth of the 
map region, while the "corner" filter transfer functions represent one-ninth of the map region located in the lower right corner. See the 
caption to Figurejsjand Section [6^ for a discussion of the o bvio us features in these transfer functions. The rotation of the features between 
center and corner transfer functions is discussed in Section 16.21 



tain an estimate of the expected depth of the maps. A 
simple calculation using the numbers for mapping speed, 
field size, and integration time to determine depth would 
yield expected depths of roughly 11 and 29 /iK-arcmin 
at 150 and 220 GHz, respectively. These numbers are 
significantly lower than the 14.8 and 36.1 /iK-arcmin cal- 
culated directly from the noise PSDs. However, several 
factors reduce the time spent on the uniform-coverage re- 
gion for which the field size of 95 deg^ is reported. First, 
with a one-degree-wide detector array we must map an 



11-by-ll-degree region to end up with a 10-by-lO-degree 
region which is seen by every bolometer. This results in 
an efficiency loss of roughly 20%. Second, the constant- 
velocity, bolometer, and scan cuts described in Section 



5.1 result in an approximately 15% loss. When these cuts 
and inefficiencies are taken into account, the expected 
depths are 13.4 and 34.7 /iK-arcmin at 150 and 220 GHz, 
within 10% of the observed map depth. All map depth 
estimates assume that the noise is white, which is known 
to be only an approximation even for the I range used 
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(a) Noise PSD for the 150 GHz map 
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(b) Noise PSD for the 220 GHz map 



Fig. 10. — Estimated two-dimensional noise PSDs for the 150 GHz and 220 GHz maps processed with point-source masking and the 
Sanson-Flamsteed projection. Note that these are averaged over the whole map area. The high-noise area near ky = kx = is due to 
atmospheric noise, the uniformly increased noise at low kx is due to low-temporal-frequency noise uncorrelated between detectors, and 
the areas of excess noise at ky ~ 15000, kx — are due to upmixed atmospheric power resulting from different sampling and filtering of 
atmospheric noise in diff'erent detectors' time-ordered data. As in the filter transfer functions (see Figures Js] and [oj, the dark stripe along 
kx = is due to the effective scan-direction high-pass filter from the polynomial and Fourier-mode subtraction, anathe hexagonally spaced 
dark spots are due to the sensitivity of the wedge mean subtraction to modes with wavelengths and angles corresponding to the detector 
array configuration. 

for these estimates. 



6.4. Comparison to Past SPT Analyses 

The maps and data products presented in this re- 
lease are very similar, but not identical, to maps of 
the ra5h30dec-55 field used by the SPT team in 
previously published scientific results. These pub- 
lished results include catalogs of SZ-selected galaxy clus; 
ters (Staniszewski et al. 2009; Vanderlinde et al. 2010 



Williams on et al. 201 ip and emissive sources ( [Vieira 



et al.n2010| , an cT measurements of the CMB angular 
power spectrum (Lueker et al.|2010 Shirokoff et al.|20ir 



|Keisler et al. 201ip . Minor variations in data selection, 
nltering, map projection, and approaches to defining 
transfer functions and noise PSDs are expected to cause 
small changes in catalogs of objects extracted from these 
maps and to estimates of the CMB power spectrum. We 
have reproduced two past catalogs and two past power 
spectrum analyses using this map release to verify that 
the differences are negligibly small. 

Using the 150 GHz map from this release that was pro- 
duced with the Sanson-Flamsteed projection and point 
source m asking, we have r e produ ced the analysis re 
ported in I Vanderlinde et al. 



clusters presented in Vander 



3procLuc( 

poTol). Ah ra5h30dec-55 
inHeet al.f2Q10 (10 out of 21 



total clusters in that work, which also included data from 
the ra23h30dec-55 field) are clearly detected in our 
analysis with the data products from this release. The 
signal-to-noise values for these 10 cluster s in our anal- 
ysis agree with the values quoted in Van derlinde et ah] 
(I2OIO) to within 2% in the mean, with an rms scatter of 
4%. This check was performed using the estimated PSD 



presented in this release, but with the analytic transfe r 
function estimate described in [Vanderlinde et ah (2010). 
A cross-check using the simulated transfer functions and 
full two-dimensional beam functions presented in this re- 
lease produces identical results at the sub-percent level. 

Using both the 150 GHz and 220 GHz maps from 
this release that were produced with the oblique Lam- 
bert equal-area azimuthal projection and no point-source 
masking, we have reproduced the analysis reported in 



Vieira et al 
noise PSDs, 
this release. 



rans 
The 



(|2010|). This check was performed using the 
'er functions, and be ams presented in 
(|2010[) catalog is faith- 



Vieira et al. 



fully reproduced using the data pro duct s in this release 
with every 5(J source from Vieira et al.| (|2010) detected 
at > 4.4(7 in our analysis and me an signal-to-noise within 
0.1% of the [Vieira et al.| ( |2010D values at 150 GHz and 
1.0% at 220 GHz. The recove red fluxes of the sources in 
this analysis compared to the Vieir a et al.| (^2010) fluxes 
are 3.1% lower at 150 GHz and 7.5% lower at 220 GHz. 
These differences are due to improved beam and calibra- 
tion estimates. If we repeat the analysis using the maps, 
fllter transfer functions, and noise PSDs p resented here 
and th e beam and calibration used in the 'Viei ra et al.l 
(2010) analysis, we reproduce the 150 and 220 GHz flux 
numbers to better than 1%. 

Finally, we construct a CMB power spectrum estimate 
from the 150 GHz data and comp are it to the power spec- 
tru m of this fleld as me asured in 'Shirokoff et al. (2011) 
(2011). In each comparison, we use the 
IS hirokoff et al.[ ( |2QlT[ ) or [Keisler 
We construct the power spec- 



and Keisler et al. 

calibration used In the 



et al. (2011) analyses. 
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(a) Center noise PSD for the 150 GHz map (b) Corner noise PSD for the 150 GHz map 
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(c) Center noise PSD for the 220 GHz map (d) Corner noise PSD for the 220 GHz map 



Fig. 11. — Estimated two-dimensional noise PSDs for selected sub-regions of the 150 GHz and 220 GHz maps processed with no 
point-source masking and the oblique Lambert equal-area azimuthal projection. The "center" PSDs represent the innermost one-ninth of 
the map region , while the "corner" PSDs represent one-ninth of the map region located in the lower right corner. See the caption to Fig ure 
[lOland Section |6.3.1| for a discussion of the obvio us features in these PSDs. The faint lines at constant kx (not apparent in the Figure [Tol 
FSDs) are due to readout lines (see Section [6.3.11 for details). The rotation of the features between center and corner PSDs is discussed m 
Section l6.3.1l 
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trum by Fourier transforming the point-source-masked 
150 GHz map ti mes an adapted version of the apodiza- 



tion mask used in 



Shirokoff et al. (2011) or Keisler et al 



( 2011 ) , correcting for the beam and transfer function, 



squaring, azimuthally averaging in t space, and subtract 
ing off a noise bias calculated by azimuthally averag- 
ing the squared noise PSD (also corrected for beam and 
transfer function). Fourier modes with significantly el- 
evated noise or where the transfer function is near zero 
were excluded from the signal and noise averages. 

We find that this simply calculated power spectru m of 
the ra5h3Qdec-55 map presented here, using the |Shi- 



rokoff et al. (2011) calibration and mas ks, agrees with 
the po wer spectrum of this field from the [Shirokoff et al. 
(2011) analysis to within 3% in power (1.5% in temper- 



ature) . This small difference is attributable to mode- 
mixing and window-function effects that are not taken 
into account in the simple analysis performed here. We 
find that the power spectrum of the ra5h3Qdec-55 map 
presented here, calculated using the Keisler et al. (2011) 



calibration and ma sks, agrees with the power spectrum of 



this field from the [Keisler et al. (2011) analysis to within 
1% in power (0.5/0 in temperature). 

7. CONCLUSIONS AND AVAILABLE DATA PRODUCTS 

We have presented the first publicly released maps 
from the South Pole Telescope, of a 95 deg^ field observed 
during the 2008 season in two frequency bands (150 and 
220 GHz). We have described the observations, data 
selection, filtering, and mapmaking approaches used to 
create these maps. In addition, we have characterized 
the instrument bands and beams, the filter transfer func- 
tions, and the noise properties of the maps. The maps 
and auxiliary data products documented in this paper 
are available online for download and use by the broader 
community. The full set of data products, available at 
| http://pole.u chicago.edu / public / data/ maps / ra5h30dec- 1 
|55" and from the NASA Legacy Archive for Microwave 
Background Data Analysis server, includes: 

• Maps: Two versions of both the 150 GHz and 

220 GHz maps of the ra5h30dec-55 field are 
given. One version uses the Sanson-Flamsteed 
projection and point-source masking in the filter- 
ing, and the other version uses the oblique Lam- 
bert equal-area azimuthal projection with no point- 
source masking. 

• Beams: Two-dimensional beam functions for 
150 GHz and 220 GHz are provided, as well as 
one- dimensional beam averages in real and Fourier- 
space representations and an estimate of beam un- 
certainties. 



• Bands: Measured band-pass functions for 150 GHz 
and 220 GHz are provided. 



• Filter Transfer Functions: Two sets of transfer 
functions are provided, corresponding to the two 
choices of projection and filtering. For the oblique 
Lambert equal-area azimuthal projection, transfer 
functions estimated for 9 sub-regions of the map 
are given. 

• Noise PSDs: Two sets of noise PSDs are pro- 
vided, corresponding to the two choices of projec- 
tion and filtering. For the oblique Lambert equal- 
area azimuthal projection, PSDs estimated for 9 
sub-regions of the map are given. 



An example of a calculation using these data products 
is described in Appendix [A] 

We have checked that these data products, when used 
in analyses similar to published analyses of this field's 
data, faithfully reproduce the published results. The 
data release and this accompanying paper are the first 
step toward an eventual release of data from the full 
2500 deg^ three-band SPT-SZ survey. 
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APPENDIX 



A. EXAMPLE CALCULATION 

As a simple introduction to the use of the data products presented here, and as a means for users to check that they 
are interpreting the data as intended, we present an example of a typical calculation, namely the estimation of the 
fiux of a point source in the map. 

The most straightforward way to estimate the fiux of a point source with a known position is to read off the map 
value at the pixel that corresponds to that position, and to convert that value in K-CMB to a value in MJy/sr using 
the conversion factors in Table [l] To convert that to a point source ffux in Jy, one needs an effective solid angle. It 
is tempting to simply use the solid angle of the beam; however, that implicitly assumes that signals in the map have 
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undergone no filtering beyond beam smoothing, whi ch is not necessari ly the case. In the case of a filtered map, the 
effective solid angle to use in this calculation is (e.g., Vieira et al.||2QlQ ) 



47r2 



(Al) 



where B{kx^ ky) is the two-dimensional Fourier- space beam function and Fikx^ky) is the two-dimensional Fourier-space 
filter transfer function. 

This calculation will give an unbiased estimate of the source flux, but the estimate will be noisier than it has to 
be, because all spatial modes are treated equally, regardless of noise. If the map is filtered further to downweight 
noisy modes, the signal-to-noise on objects with a known angular profile (such as point source s) can be improved. The 
Fourier-space optimal filter for sources of a known shape is (e.g., Haehnelt fc Tegmark||1996 ) 



( ^ ky^N ( kx 1 ky ) 

\/ ^'^{kx-) ky)^ {kx-) ky) 



(A2) 



where r is the (Fourier-space, beam-and-filtering-convolved) source profile, and N is the (Fourier- space) noise covari- 
ance matrix, which we assume to be diagonal. The noise covariance can include contributions from unwanted sky 
signals (convolved with the beam and filtering) as well as instrumental and atmospheric noise. For point sources, the 
source profile is entirely a function of the beam and filtering: 



rps{kx, ky) = B{kx, ky) F{kx, ky) 



(A3) 



To extract source fluxes from the optimally filtered map, the values at the pixel location should be multiplied by the 
conversion factor in Table [l] and an effecti ve so lid angle that now also includes the optimal filter tp; however, if the 
filter is properly normalized as in Equation |A2[ then the effective solid angle will not change. 

We can now estimate the flux of a known source in the SPT map using the raw data and an optimally filtered map. 
The radio source PKS 0549-575 (also known as SPT-S 055009-5732.3) is located at R.A. 87.5399°, decl. -57.5401°, 
which corresponds to pixel location [937,928] (zero-based indexing) in the oblique Lambert equal-area azimuthal 
projection maps. The values in the point-source-unmasked 150 and 220 GHz maps at this pixel location are 10.084 x 
10~^ and 8.508 x 10~^ K-CMB. Multiplying by the radio-source conversion factors in Table [l] yields values of 3.996 and 
4.056 MJy/sr. Multiplying the Fourier-space beam by the filter transfer function appropriate to the source location in 
each band and integrating over all k values yields effective solid angles of 1.523 x 10~^ and 1.220 x 10~^ sr, leading to 
source flux estimates of 0.6086 and 0.4949 Jy. We can estimate the statistical uncertainty on these values by measuring 
the pixel variance in the neighborhood of this source and converting that number to Jy. This procedure yields rms 
uncertainties of 0.0043 and 0.0081 Jy at 150 and 220 GHz. 

If we assume that the only signals in the map are CMB fluctuations and point sources, we can build an optimal 
filter using the source profile defined in Equation A3 and a noise power spectrum that is the sum of the (squared) 
instrument-plus- atmosphere PSD and the CMB power spectrum. Because we have filtered the CMB, the version of 

shows the 



Al 



the power spectrum that goes into the optimal filter must be multiplied by [B{kx^ky)F{kx^ky)]'^ . Figure [ 
azimuthally av eraged optima l filter at 150 and 220 GHz using the data products presented here and the (jMB power 
spectrum from Larson et al. (2011). After convolving the maps with these filters, the map values at the location of 
PKS 0549-75 are x l0^"aiid 8.761 x 10"^ K-CMB at 150 and 220 GHz. As expected, the effective solid angle 

is the same as before application of the optimal filter, and the source flux estimates obtained by multiplying these 
optimally filtered map values by the radio-source conversion factors in Table [l] and the solid angles calculated above 
are 0.6194 and 0.5096 Jy. Estimating the flux uncertainties as we did with the unfiltered map yields rms uncertainties 
of 0.0016 and 0.0032 Jy at 150 and 220 GHz. For a source this bright, these statistical uncertainties are dwarfed by 
the calibration and beam uncertainties. 
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